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SUMMARY 


This work is concerned with large deflections of a shallow elastic conical 
membrane fixed at the outer edge and loaded by either uniform or hydrostatic 
pressure. The governing equations were solved by the method of matched asymp- 
totic expansions and by a finite- difference method. Agreement between the two 
methods was excellent for the small values of the perturbation parameter. 


INTRODUCTION 


This paper is concerned with the moderately large axisymmetric deformation 
of a shallow, elastic conical membrane. The purpose of this work is to further 
investigate the application of the method of matched asymptotic expansions (see 
Van Dyke, reference 1) to the solution of membrane-shell problems involving 
large deflections. The success of this method is based on the fact that for 
small loads the linear membrane solution is a good approximation to the actual 
solution everywhere except in the immediate vicinity of boundaries. In these 
regions thin boundary layers exist where the variables undergo rapid changes to 
accommodate themselves to the boundary conditions that cannot be satisfied by 
the linear membrane solution. In the method of matched asymptotic expansions 
separate perturbation expansions ar-e found in the interior and boundary- layer 
regions and matched in an appropriate way to insure that they join smoothly. 

Bromberg and Stoker (ref. 2) initiated this type of analysis of membrane 
shells when they found one term of both the interior and boundary- layer expan- 
sions for a uniformly pressurized shallow spherical shell. The next two terms 
in the interior and boundary-layer expansions were found by Smith, Peddieson, 
and Chung (ref. 3) and used by them to investigate the accuracy of finite- 
difference solutions of the same problem. One term of the interior and boundary- 
layer expansions for deep membranes of arbitrary shape has been given by 
Rossettos (ref. 4). This work generalizes the results given in the references 
listed in reference 4. 

In the present paper three terms of the interior and boundary- layer expan- 
sions are found for the ease of a shallow conical membrane loaded by either 
uniform or hydrostatic pressure. It is found that complications arise which do 
not appear in the solution of the corresponding sphere problem. The solution 
method is modified somewhat to account for this. Numerical results are pre- 
sented to illustrate some of the interesting features of the solution. 


575 



GOVERNING EQUATIONS 


Consider a shallow conical membrane (opening upward) with base radius a, 
thickness h, and initial angle <{io with the horizontal made of a linearly elastic 
material with modulus of elasticity E and Poisson's ratio v. The equations gov- 
erning moderately large axisymmetric deflections of such a structure can be 
obtained from the work of Reissner (ref. 5). The resulting equations are (in 
dimensionless form) 

^ + ip'/r - + (l+e^6/2)(B/r) = 0 

(l+E^g);!, = rV • 

N^ = ij;/r, Nq = li;’ 

u = r;/)' - vi/j, w' = B (1) 

where ar is the radial coordinate, VQaiii/<t)o is a stress function (Vq being a 
characteristic vertical force resultant), VqV is the vertical force resultant, 
VoNp/<j>o is the radial stress resultant, VoNq/^q is the transverse stress resul- 
tant, e is a load parameter, atfios^u is the horizontal displacement, a(j)GS^w is 
the vertical displacement, is the rotation, and a prime denotes differen- 

tiation with respect to r. 

In the present paper a uniform pressure po and a hydrostatic loading 
Yo4>o3-(l”r) are considered. It can be shown by considering the vertical equi- 
librium of the membrane centered on the vertex and having radius r that 

V = r/2 - jr2/3 (2) 

where j a 0 for the uniform pressure and j = 1 for the hydrostatic pressure . 

The characteristic vertical force resultant is given by 


V 

o 


{ 


p a 
o o 
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j = 0 

j = 1 


The load parameter e is defined to be 

e = (V /Eh<f>^)^ 
o o 


(3) 

(4) 


In the present work it is desired to solve equations (1) subject to the 
boundary conditions 


u(l) = w(l) = 0 


(5) 


Special attention will be given to situations where e << 1. 
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STRAIGHTFORWARD SOLUTION 


To begin the solution process a straightforward perturbation solution to 
equations (1) is sought for e « 1. To do this it is convenient to rearrange 
equations (la) and (lb) to yield 

e^(i/;" + i|/V/r - ip/r^) - (1 - (rV/i|()^)/(2r) = 0 

S = (rV/iJ) - l)/e^ (6) 

A straightforward perturbation solution for e << 1 has the form 


^SO 


si 


s2 


(7) 


where the subscript s indicates the straightforward solution. Substituting 
equation (7) into equation (6a), expanding for e << 1, setting the coefficient 
of each power of e equal to zero in the usual way, and solving the resulting 
algebraic equations yields 

jp ^ (r^/2 - jr^/3) + e^(r^/2 - jr^/3)(3r/2 
s 

- 8jr^/3) + e‘*^(r^/2 - jr^/3)(75r^/8 

- j(79r3/2 - 32r^)) + . . . (8) 

From equations (1) and (6b) it can then be shown that 

<v (r/2 - jrV3) + e^(r/2 - jrV3)(3r/2 

- djV‘^/3) + e‘*(r/2 - jrV3)(75rVs 

- j(79rV2 - 32r‘^)) + . . . 

N % (r - jr^) + e2(9rV4 - j(22rV3 - 40rV9)) 

0S 

+ e^(75r ^/4 - j(915rV8 - 175r^ 

+ 224rV3)) + . . . 

8g '\--(3r/2 - 8jr^3) - e^(57r^/8 

- j(63r^/2 - 224rV9)) + . . . 

% (1 - v/2)r2 - j(l - v/3)r3 + e^(3(3 - v)r^4 

s 

- j(ll(4 - v)rVe - 8(5 - ^^)rV9)) 

+ e'+(75(4 - v)rVl6 - j(183(5 - ^)rVs 

- 175(6 - v)rVe + 32(7 - v)rVs)) + . . . 
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w -3(r^ - l)/4 + j8(r^ - l)/9 - e^(19(r^ - l)/8 
s 

- j(63(r‘* - l)/8 - 224(r‘^ - l)/45)) + . . . 


( 9 ) 


where equation (5b) has been used to determine the constants of integration in 
equation (9c). By comparison with the results given in Kraus (ref. 6) it can 
be seen that the first term in each series expansion is the linear membrane 
solution. It should also be noted that the first terms in equations (9d) and 
(9e) are due to the second term in equation (8). Thus to obtain 8 and w to 
O(e^) it is necessary to find to 0(e‘^). The boundary condition represented 
by equation (5a) cannot be satisfied by equation (9d). Thus a boundary- layer 
expansion is needed in the vicinity of r = 1. 


BOUNDARY-LAYER SOLUTION 


There are several ways to carry out the b oundary- layer analysis in this 
problem. One is to work in terms of the original stress function If this 

is done the differential equation foi' the first boundary- layer approximation 
turns out to be nonlinear . Bromberg and Stoker (ref. 2) discovered that a 
linear equation could be obtained in the first approximation for a spherical 
membrane by a method which is equivalent to working with a dependent variable 
which is the difference between the actual and the linear stress functions . 

This was tried in the present problem but matching difficulties were encountered. 
These were due to the fact that equations (8) and (9) do not terminate with one 
term for the cone as the corresponding straightforward expansions do for a 
sphere. It was, therefore , decided to use the difference between the actual 
stress function and the straightforward stress function as the dependent vari- 
able. This guarantees that the outer expansion for this dependent variable will 


be zero. Thus it is necessary to find only the inner expansion. 

Substituting 

>1' = ’I's ’^b (10) 

(where the subscript b denotes the boundary- layer solution) into equation 
(6), defining the boundary- layer variables F and 5 by the equations 

li'jj - eF, r = 1 - eg, (11) 

expanding F as 

F 'V F^ + eF^ + e^F^ + . . . , (12) 

and carrying out the usual perturbation analysis yields 

Fjj - s2Pjj = 0 (13) 

and two other equations governing F^ and where 

S = (6/(3 - 2j))^ (14) 
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In equation (13) ( ) = d( )/d^. Define Ngj^, Uj^, and Wj^ by the follow- 

ing equations 


N = N + eN , , 
r i>s rb’ 


“e = “es ”eb 


B = 3„ + 6t,/£» u = u + 


s ^b 
Now expand as follows 




w = w„ + w, 
s b 


(15) 


'S, "" *bo ®*bi + * • • • 

where A]-, is any one of the bo\mdary- layer variables. Substituting equations 
(10), (11), (12), (15), and (16) into equations (l) , ejq>anding for e « 1, and 
equating the coefficients of like powers of e to zero one obtains 


N • 

rbo 

= V 

“ebo = ®bo = 


11 

O 


”ba = 

(17) 


and two similar sets of equations relating and A, ^ to F^, , and F^. A 

similar procedure applied to equation (5a) leads to boundary condition 

Fq( 0) = 1 - j - (1/2 - j/3)v (18) 

• • 

and boundary conditions for F^(0) and F2(0). 

To illustrate the solution procedure the first approximation will now be 
carried out in detail. The solution of equation (13) is easily seen to be 

F^ = Cj exp(S|) + c^ exp(-S?) (19) 

Since the outer expansion has been forced to vanish because of equation ( 10 ) 
the matching process (see Van Dyke, reference 1) is equivalent in this case to 
a statement that positive exponential terms must vanish. Thus 


c^ = 0 (20) 

Substituting equations (19) and (20) into equation (18) yields 

c^ = -(1 - j - (1/2 - j/3)v)/S (21) 

Thus 

F^ = -(1 - j -(1/2 - j/3)v)exp(-Se)/S (22) 

Siibstituting equation (22) into equations (17) one obtains 
Nrbo " "(1 “ 9 -(1/2 - j/3)v)exp(-S5)/S 
Ngbo = (1 ■ i ■ (1/2 - j/3)v)exp(-S5) 
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$bo = (1 - i - (1/2 - j/3)v)S exp(-sc) 

= -(1 - j - (1/2 - j/3)v)exp(-S?) 

= -Cl - j - (1/2 - j/3)v)(l - expC-Sg)) (23) 

The results for higher approximations are found in a similar way but the calcu- 
lations are quite lengthy. For the sake of brevity this work is omitted. 

To find the complete solution the boundary- layer expansions must be added 
to the corresponding straightforward expansions. The first approximations to 
these expressions are 

- r^/2 - jr^/a - e(l - j - (1/2 - j/3)v)exp(-S(l - r)/e)/S 

= r/2 - jr2/3 - e(l - j - (1/2 - j/3)v)exp(-S(l - r)/e)/S 

Nqq = r - jr^ + (1 - j - (1/2 - j/3)v)exp(-S(l - r)/e) 

3q = (1 - j - (1/2 - j/3)v)S exp(-S(l - r)/e)/e 
Uq = (1 - v/2)r2 - j(l - v/3)r^ - (1 - j 

- (1/2 - j/3)v)exp(-S(l - r)/e) 

w = 3(1 - r2)/4 - 8j(l - r3)/9 - (1 - j 
0 

- (1/2 - j/3)v)(l - exp(-S(l - r)/e)) (24) 

In writing equations (24) the bcundary- layer solution was treated as the funda- 
mental expansion. All terras in the straightforward expansion with magnitude 
equal to or greater than the first term in the boundary- layer expansion were 
added to this terra to form the first approximation. The same method was used 
to obtain the second and third approximations. 


RESULTS AND DISCUSSION 


Numerical results were computed for the first, second, and third approxi- 
mations to the variables ip, Np, Nq, 3, u, and w. These calculations were made 
for a variety of values of the load parameter e and Poisson’s ratio v. To 
evaluate the accuracy of the perturbation method, selected cases were compared 
with numerical solutions to equation (6a) obtained by the finite -difference 
method discussed by Smith, Peddieson, and Chung (ref. 3). It was found that 
the third approximation to the peidrurbation solution agreed with the finite- 
difference results up to e = 0.1. It should be pointed out that for small 
values of e,. the numerical method is difficult to apply because a variable 
step size must be used near the edge and the optimum arrangement of step sizes 
can only be approached by trial and error. The explicit formulas obtained in 
the present work are much easier to use for e « 1. 
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To illustrate the behavior of the solution some of the computed results 
are shown in figures 1 --• 4 . For the sake of brevity, data are presented 
for only the radial stress resultant Np, the transverse stress resultant N9, 
and the vertical deflection w. The solid lines represent the three-term per- 
turbation solution while the dashed lines represent the linear membrane solu- 
tion. The linear membrane solution is shown only when it differs significantly 
from the perturbation solution. 

Figures 1 and 2 present results for uniform pressurization (j =0). 
Figure 1 shows that thin boundary layers exist for N9 and w for e = 0.01 
while Np does not exhibit boundary-layer behavior. As e incr’eases the bound- 
ary layers become wider for all variables. This is illustrated by figure 2. 
Figures 3 and 4 contain results for hydrostatic loading (j = 1 ). The para- 
metric trends illustrated by these results are identical to those discussed 
above but the behavior of the solution variables is more complicated. These 
results illustrate the utility of the perturbation method. Complicated func- 
tions of this type can be represented numerically only if extreme care is used. 

Results were also computed for several other values of v. It was found 
that the qualitative behavior of the solution is not significantly influenced 
by this parameter. 


CONCLUSION 


In this paper, the rotationally symmetric moderately large deformation of 
a linearly elastic shallow conical membrane siibjected to either uniform or 
hydrostatic pressure was investigated. A single differential equation having 
a stress function as dependent variable was solved by the method of matched 
asymptotic expansions. The accuracy of the solution was verified by compari- 
son with a finite-difference numerical solution of the governing equation for 
the stress function. Selected results were presented graphically to illustrate 
interesting features of the solutions . 
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Figure 3*“ Stress resultants and axial Figure 4.— Stress resultants and axial 

displacement for hydrostatic loading, displacement for hydrostatic loading 






